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We present a new mechanism for generating primordial statistical anisotropy of curvature per- 
turbations. We introduce a vector field which has a non-minimal kinetic term and couples with a 
waterfall field in hybrid inflation model. In such a system, the vector field gives fluctuations of the 
end of inflation and hence induces a subcomponent of curvature perturbations. Since the vector 
has a preferred direction, the statistical anisotropy could appear in the fluctuations. We present 
the explicit formula for the statistical anisotropy in the primordial power spectrum and the bispec- 
trum of curvature perturbations. Interestingly, there is the possibility that the statistical anisotropy 
does not appear in the power spectrum but does appear in the bispectrum. We also flnd that the 
, statistical anisotropy provides the shape dependence to the bispectrum. 

(N 

I. INTRODUCTION 

, Current cosmological observations such as cosmic microwave background (CMB) observations tell us that our 
C^' universe is likely to be homogeneous and isotropic. However, there is the possibility that the effect of rotational 
, invariance violation has been already detected in the CMB ternperature anisotropies, as an apparent alignment of the 
'. CMB multiples on very large scales, so-called " Axis of Evil" [E i, i, i, i, i, 0, i, i, E [H H [H, [3, 111, d, [13, H 
' [l^, [13, Hll, [22] ■ If it is true, the statistical anisotropy might stem from the physics in the early universe. Hence, the 
Ci . possibility of statistical anisotropy deserves further investigation. 

On the theoretical side, statistically homogeneous, isotropic, and Gaussian density fluctuations had been supposed 
, to be robust prediction of the inflationary scenario. However, the possibility of the statistical non-Gaussianity has 
J> ■ been studied extensively. It turns out that it is not difficult to make a non-Gaussianity if we incorporate the multi 
' scalar fields into the system. In fact, there are many scalar fields in any realistic particle physics model. Hence, it is 
^^O natural to test the non-Gaussianity seriously by making use of observational data. For the same reason, in this paper, 
■ we would like to argue that the statistical anisotropy should be checked observationally. The point is that there are 
many vector fields in addition to the scalar fields in the particle physics model. As we will see, the vector field could 
violate the rotational invariance at the end of inflation if it couples with a waterfall fleld in a hybrid inflation model. 
^5 ■ That means the statistical anisotropy could be generated due to the ubiquitous vector fields. Of course, the curvature 
perturbations induced by the vector field should be subdominant in the total curvature perturbations in order to 
reconcile the scenario with the current data. Nevertheless, the effect of the vector field should not be overlooked 
under the current precision cosmology. 

The paper is organized as follows. In section [TTl we show a new mechanism generating the statistical anisotropy 
at the end of inflation. This model is based on Ref [H, Hi, HE [H, [13, HI, [H , in which the authors consider the 



..... . 

modified hybrid inflation adding a light scalar field other than inflaton in order to generate the large non-Gaussianity 

of primordial curvature perturbations. We put a massless vector field instead of a light scalar field. The vector 
field violates the rotational invariance and as a consequence generates the statistical anisotropy in the curvature 
perturbations. In section lllli we present the calculation of vector fiuctuations. It allows us to explicitly show 
the effects of violation of rotational invariance on the power spectrum and the bispectrum of primordial curvature 
perturbations. We also provide a simple example as an illustration. We find that the statistical anisotropy in the 
bispectrum could be more significant than that in the power spectrum. In section ITVl we give a brief summary. 
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II. BASIC IDEA 

In this section, we give a basic idea of the new mechanism for generating the statistical anisotropy in primordial 
curvature perturbations. The explicit calculation will be given in the next section. 

In the original hybrid inflation model, when the effective mass squared of a waterfall field becomes negative, i.e., 
the waterfall field becomes tachyonic, the tachyonic instability makes the waterfall field roll down toward the true 
vacuum state and the inflation suddenly ends. The critical value 4>e, which is the value of the inflaton at the end 
of inflation, should be determined by the model parameter and hence 0e does not fluctuate. However, if some other 
fields couple with the waterfall field, the situation would be different. 

Let us consider a system constructed by two scalar fields, say, an inflaton (f) and a waterfall fleld %, and a vector 
field Afj,{ii = 0, 1, 2, 3) which couples with the waterfall field. The action can be written as 

-\ J d\^gg^^'^gP^f{^)F^,F,, , (1) 

where F^^ = d^ A^ — duA^ is the field strength of the vector field, V{(f), A^) is the potential of fields, and an arbitrary 
function /(</)) represents gauge couphng. We have to choose V{(j), x-> Ap) so as not to violate the gauge invariance. For 
the background universe to be isotropic, we need the condition that the energy of the inflaton should be dominant in 
the total energy of the universe during inflation era. Under this condition, we have a small expectation value of the 
vector field. Here, we do not specify the form of potential but we assume that the vector field is massless and have a 
small expectation value compared with the inflaton. Hence, we neglect the terms which come from the coupling with 
the vector field in the background equations of motion for the scalar fields and we treat the vector field perturbatively. 

A. formalism 

First, we briefly review the bN formalism [sO, HH, [12, HI] (or separate universe approach [13]), which is well-known 
as the most powerful tool to calculate the evolution of curvature perturbations on super-horizon scales. Based on the 
formalism, the curvature perturbations on super-horizon scales, where spatial gradients and anisotropic shear are 
small, can be written by 

C(ti=^)~5iV(iF,4) , (2) 

where SN{tp^ i») represents the perturbation of the e-folding number N measured in homogeneous background universe 
between the initial flat hypersurface at i = and the final uniform energy density hypersurface at t — tp. More 
precisely, the e-folding number is given by iV = J^*^ Hdt , with the Hubble parameter H . When we take the final 
time tp to be a certain time after the complete convergence of the background trajectories in phase space, C, becomes 
constant. Let the initial time be some time during the slow-roll inflation driven by the inflaton. Soon after the time 
that the scale of interest crosses the horizon, we can rewrite Eq. ([2]) as a Taylor expansion up to the second order 

atp) = NiMi + \Ni.jM'Mi , (3) 

where Ni^ = dN{tp,t^) / d4>^ {t,:) and A^/j* = d'^N{tp,t^)/d(p^ {t,,)d(j)'^ (t^,). Here, we consider that the background 
trajectories are parameterized by the scalar field (f)^ . Hence, the perturbations of the scalar field Scp^ give rise to the 
difference between the evolutions of the adjacent background homogeneous trajectories. 

B. Curvature perturbation generated at the end of inflation 

As we mentioned in the previous subsection, the curvature perturbations on super-horizon scales are given by the 
perturbations of e-folding number SN . In the standard single scalar inflation or hybrid inflation, the inflation ends 
when the inflaton cj) reaches a critical value (jje- On the other hand, in the multi-component inflation, the critical value 
(j)e may fluctuate due to a light fleld other than the inflaton (p. Hence, in such a situation, there is a possibility of 
generating the curvature perturbations through the fluctuation of 0e . Based on the models introduced by Ref. [23[ . , 
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in this paper, we introduce a massless vector field as another hght field. Using the 5N formula, the curvature 
perturbation generated at the end of inflation {t = t^) can be expressed as 



dN 1 d^N 



dN dMA) ^ 1 



2 d^l 

ON 0e {A) d^N d(t>e (A) d<j)e (A) 



d(j>e d Aid AO d(i)l dAi dA^ 



5K5Ai , (4) 



where we set A^ — 0. That is, 0e can fluctuate due to fluctuations of the vector field A*. Let us take the hypersurface 
at the end of inflation t = tf. io be that of a uniform energy density. Then, the total curvature perturbations at the 
end of inflation t = is given by 

C(^e) = Cm/ + (,end , (5) 

where 

dN ^ 1 92 iv „ 
C„, = — + . (6) 

We assume that all fields decay into radiation component through the sudden reheating just after the end of inflation. 
Under this assumption, the evolution of the universe can be described by one-parameter, namely, radiation energy 
density, after the end of inflation. Hence, the curvature perturbations on a uniform energy density hypersurface 
become constant after the end of inflation and the iso-curvature mode does not persist until present. At the leading 
order, the power spectrum of curvature perturbations is given by 

(Ck.Ck.) = Pc(ki)5^^^ (ki +k,) = NlP,{k^)5(^^ (k, +k,) +iV|^M^^M^(M:,/Ai,^) , (7) 

where A^* — dN/d(l)t., — dN/d(j)e and (i^^ki'^^k,) power spectrum of the perturbation of the vector field, 

which we will calculate in the next section. Here, we neglect the cross terms such as {6(j)^,u-^6Aly^^) and assume that 
(5(/>* has pure Gaussian distributions. To the leading order, the bispectrum is also given by 

(CkiCk.Ck3> = (2^)~3/'i3c(ki,k2,k3)<5(3) (ki+ka+kg) 

= (27r)-3/2iV2Ar,, [P^{ki)P^{k2) + 2 perms] 5^3) (ki + kz + kg) 

, dMA) d<PM) ^'^e(^) , , , , 

dA^ dAo dA^ (-^^^k/^^k.-^^^ks) 



^dMA)dMA)f^ d^eiA) ^ NeedMA)dMA) 



dA^ dAi KNedA'^^dA^-^ dA^-^ dA'^^ 

X '(SAl^JAi^^ {SA'^^SA'^)^J + 2 perms 



(8) 



where TV** = 92iV/90*9(/), , = d'^N/d(f)ed(f>e and ★ denotes the convolution. 

From these expressions, we find that one can obtain the curvature perturbations which have the direction- 
dependence due to the vector field perturbations. 

III. PRIMORDIAL STATISTICAL ANISOTROPY 

In this section, we first quantize the vector field with a non-minimal kinetic term and obtain the power spectrum of 
the vector fluctuations. Next, we derive the explicit formula for the power spectrum and the bispectrum of curvature 
perturbations. Then, we present a simple example as an illustration. 

A. The spectrum of vector perturbations 



Here, we give the spectrum of vector perturbations, following Rcf. [35], which appear in the expression of the 
spectrum of curvature perturbations ([7]) and ^ . 
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Following Ref. [3a|, we adopt the Coulomb gauge where 5Aq{x'^) = and djA^{x'^) = 0. Then, we obtain the 
evolution equation for the linear perturbation as 



A'^ - J-A^ - a^^J^^A^ = 



(9) 



where Ai = f6Ai, a prime represents the derivative with respect to a conformal time, r = J dt/a, and we neglected 
the potential term. Here, we assumed the gauge invariance for the vector field. Hence, there should be a derivative 
coupling between the vector field and the waterfall field. However, this term is not relevant to our analysis. 
Now, we can expand the vector field in terms of annihilation and creation operators, dAk and a\y^, as 



J2 e»A (k) Vk (T)aAke*-" + vl iT)al^^e 



where k is the comoving wave number and the transverse polarization vectors eix defined by 



(10) 



(11) 



are introduced for a consistent quantization in the Coulomb gauge. With the change of the variable tpk = oiVk, the 
equation of motion for the mode function ijjk is given by 



f 



Using this variable, we obtain the power spectrum of vector perturbations as 

(M.(k)M,(k')) 



2 f 2 \ 



^(3)(k + k') 



(12) 



(13) 



where ki = Sijk^ . In the short wavelength limit /c^ — > cxd, from Eq. (jl2p . we see one can take the Bunch-Davies vacuum 



2k 



(14) 



Let us consider the case in which the vector perturbations SAi stay constant in time on super-horizon scales. From 
Eq. (I13p . we have to obtain jV'fcl oc a(T)/(r) in order to achieve such situation. In the slow-roll limit (e = —H/H^ 0), 
we have 



(15) 



Hence, we assume that dr = da/{o?H) and assuming that the form of the gauge coupling function f oc a", we have 



2 a{a + l) 



With the Bunch-Davies vacuum p^ . this equation can be solved easily as 



(16) 



(17) 



where hI^\x) is the Hankel function of the first kind. If we take a = 1, in the long wavelength limit — fcr 0, we 
obtain the asymptotic solution for as 



^fc(a = 1) 



{^ikr) ^ 



2k 



(l + 0((-fcr)2)) 



(18) 
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Hence, we achieve (SAiSAj) 
spectrum 



constant in time on super- horizon scales {k/aH — > 0). Thus, we have the vector power 



2k^fS 



_^]^(3)(k + k') 



fc2 



(5(3) (k + k') 



(19) 



where we defined P</,(fc) = H^/2k^. Strictly speaking, in order to evaluate the bispectrum we need to calculate 
(SAlj^^SAl^^^SAl^^^). However, as we saw in this subsection, initially the vector perturbation in Coulomb gauge can be 
considered as the Gaussian perturbation as well as the perturbation of scalar field in the standard slow-roll infiation 
model. In the standard slow-roll inflation model, as is well known, three point correlation function of the perturbation 
of scalar field is suppressed by the slow-roll parameters [36.] , that is, we can consider that the perturbations of scalar 
field at the horizon crossing time t = t^, almost obey Gaussian statistics. In this model, the vector field can be 
considered to be slowly rolling, so the deviation from Gaussian statistics of the vector perturbations 6A^ at < = 
would be too small to be detected in the future experiments, e.x., Planck [s^jUll. Hence, we assume, here, that the 
non-Gaussian feature of the vector field perturbation at i = is negligible. 



B. Rotation Violating Power spectrum and Bispectrum 



Now, we are in a position to write down the formula for the power spectrum and the bispectrum of curvature 
perturbations. 

Using Eq. (|19p . we can rewrite the expression for the power spectrum of curvature perturbations in the modified 
hybrid infiation model ([7]) as 



(CkXk.) =Pc(ki)<5(3) (ki+k2) 



Nf 



NA^ k\ki 



kl 



P4ki)6'^^^ (ki+k2) 



(20) 



where we have defined qi = d(j)e/dA^, qij = cP(f>e/dA'^dA^ . In order to compare this expression with that of Ref. ]la^ . 
can be rewritten as 



Pc(k) = PrW [l + .9/3(q-k) 
where we used the notation q = q/ |q|, k = k/ |k| and the isotropic part 

is separated. The coefficient of the anisotropic part reads 



913 



1 + /3' " \NJ, 
where we defined N^, — l/y/2e^. Here, e = V? /2V'^ is a slow-roll parameter. 



(21) 



(22) 



(23) 



The bispectrum ^ also becomes 
Pc(ki,k2,k3) = 



.rznr , \ /I N^e 

N^N^*+ [—j q,qj \^—qi^i^ + -^qiiQe, 



^2 



-xP^{ki)P^{k2) + 2 permsj 



(24) 



where we assumed that 5Ai is Gaussian. Now, we define the non-linear parameter /at^ as the bispectrum normalized 
by the isotropic part of power spectrum P^'^°{k); 



6 „ u u ^ - Pc(ki,k2,k3) 

5 P^'^°[ki)P^'^°[k2) + 2 perms 



P(ki,k2)fci 



3 



2 perms 



(25) 
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where we assumed the scale-invariant power spectrum and defined 



F(ki,k2) = (1 + /3) 



_2 ^^^^ 



■90 



2 Nee 



l-(q-ki)2 l-(q-k2 



+ 



9p 



Ne Iq^ 



(q • k2)fc^^ 



(26) 



Strictly speaking, we should have defined the non-linear parameter by using the full power spectrum instead of the 
isotropic one. However, the resultant expression is too complicated to display here. When the anisotropic part in the 
power spectrum is small, the above definition gives a good approximation to the true non-linear parameter. Indeed, we 
will discuss this interesting possibility later. Otherwise, we should go back to the original formula for the bispectrum. 

Under the slow-roll approximation, we obtain N^^./N'^ ~ Nee/N^ — 0(e). Hence, neglecting the first and second 
terms in the right hand side of the above equation, wc can further simplify the expression ()26|) as 



i^(ki,k2) 



9^ 



q'^ -iq■k^)k{^ - (q . k2)fc^^ 



(27) 



where 
as 



911^2/ 111 • We can decompose the non-linear parameter into the isotropic part and the anisotropic part 

(28) 



^/jVL(ki,k2,k3) = ^/;,l(l + ,r"(ki,k2,k3)) , 



where the isotropic part reads 



NL 



and the anisotropic part is deduced as 

q^j 



r™(ki,k2,k3) 



qm^q^ 



ct fci(q • ki) + fc^2(q • k2) + k\ki{ci ■ ki)(q • k^) 



1.3 



2 perms 



(29) 



(30) 



Here, we used the relation Ne = l/-\/2eI. Taking a look at the above formula, we notice that the statistical anisotropy 
gives the shape to the bispectrum even for the local model. 



C. Simple example 

Let us consider a simple example given by 

^(0, X, A^) = ^{x'- v'f + + Im^^' + ^h^A^A^x' , (31) 

where 0, x ^-nd are the inflaton, the waterfall field, and the light vector field, respectively. The coupling constants 
are denoted by X,g, h, the infiaton mass is given by m, and the vacuum expectation value for x is represented by v. 
The effective mass term of the vector can be accommodated in the system without violating the gauge invariance. 
Indeed, the potential (31) looks like an abelian Higgs model in unitary gauge. Namely, x could be the radial direction 
of a complex scalar field which is invariant under a local U{1) gauge transformation and (p would be a gauge singlet. 
For this potential, we can realize the hybrid infiation and obtain the evolution of inflaton field, (f), as 

^ cx a-™'/^° , (32) 

where Vq = Hence, when /((/)) have f (X a . This behavior was assumed in the previous 

subsection IIII Al in order to achieve 5Ai ~ const, on super-horizon scales. The effective mass squared of the waterfall 
field X is given by 

ml EE ~Xv^ + g^cl)^ + h^A'Ai , (33) 
where we set A'^ = 0. At the end of infiation, the equality; 



Xv^ = g^cfl -h h'A'A, , 



(34) 
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is satisfied. As is expected, the critical value, (j>e, depends on A\ Hence, we have 



j11 
Ja\ 



lAI 



j5'' + q'q' 



and 



where we used the approximation 



/3: 



1 fh^\A.\ 



(35) 



(36) 



Now, we can write down the power spectrum 

Pc(k) = Pf°(fc) 



1 



tf\A\ 



r(q-k)' 



(37) 



It should be noted that the power spectrum is scale invariant although there is a direction dependence. The magnitude 
of statistical anisotropy is determined by the parameter, g/j. 

Let us move on to the bispectrum which is more interesting. Substituting the expressions (|35p and ([551) i^to 
Eqs. (p9)) and pO)) . we have the isotropic part of the non- linear parameter 
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NL 



-Ve9l 1 + ^ 

^ \ \A\' 



(38) 



where rj = V^^/V. The anisotropic part reads 

r™(ki,k2,k3) = - ' 




2 



+ 2 perms 



l3 

"■3 



+ 2 perms 



(39) 



where the relation between wavenumber vectors holds ki + k2 + ka = 0. From the above expression, we see that 
the amplitude of the non-Gaussianity is dependent on the magnitude of statistical anisotropy, gp. For large /3 
(/3 S> 1), we can approximate ~ 1, and then, the statistical anisotropy appearing in the primordial power spectrum 
(1^^ becomes large. On the other hand, for small /3 (/3 <C 1), we can approximate ~ /3 <C 1. Naively, in this limit, 
the statistical anisotropy is small and the non- linear parameter /jvl is also too small to be detected in the future 
experiments. However, Eq. (|29p tells us that there may be a possibility of generating large non-Gaussianity even for 
the small statistical anisotropy in the power spectrum if we choose a small value for and much large value for 



while keeping gp <ti 1. From Eqs. 



and 



we also find that the anisotropic part has the same order 



/l2|A|" 

as the isotropic part in the bispectrum, in contrast to the power spectrum. Hence, it may be possible to detect the 
statistical anisotropy in the bispectrum with the future experiments. It will give us information about a new physics 
in the early universe associated with the violation of the rotational invariance. 

In this model, it is important to set the initial value of the vector field to be small. In general, the vector field A'^ 
may not be trapped near the ground state. In such a case, however, the prediction contradicts the observations, that 
is, this case predicts the anisotropic universe. Hence, it is not our universe. Of course, there may be more attractive 
scenario which does not resort to the anthropic argument. 



IV. SUMMARY 



We presented a new mechanism for generating primordial statistical anisotropy of curvature perturbations. We 
have considered a vector field coupled with a waterfall field in hybrid inflation model. We have shown that, in such a 
system, the statistical anisotropy could appear not only in the primordial power spectrum but also in the bispectrum 
of curvature perturbations. From the formula for the power spectrum and bispectrum of the curvature perturbations. 



8 



we see the magnitude of the direction-dependence is related to that of non-Gaussianity of curvature perturbations. 
The degree of the anisotropy depends on the couphng parameter between the vector field and the waterfall field 
and also the coupling between the inflaton field and the waterfall field. Intriguingly, it turned out that there is the 
possibility that the statistical anisotropy does not appear in the power spectrum but does appear in the bispectrum. 
It is also interesting to observe that the shape of the bispectrum changes due to the statistical anisotropy. It is known 
that the shape of the bispectrum is different from model to model. The statistical anisotropy may have to be taken 
into account when the data is interpreted. 

Recently, many authors discussed the statistical anisotropy of CMB temperature anisotropics, so-called "Axis of 
Evil". The resulting expressions for the spectrum of curvature perturbation, (p2|) . ((29)l and ((30)) can be compared 
with observations. If this statistical anisotropy were detected significantly, the mechanism presented in this paper 
would explain such effects without contradicting other current observations, just by tuning the coupling parameters. 

In this paper, in order to obtain the spectrum of perturbations of vector field we considered a system in which 
conformal invariance is violated based on Ref . [35| and assumed that a — 1 which represents the dependence of coupling 
function / on a scale factor a. We think that other forms of function / are also possible. It would be interesting to 
investigate those cases. We can also consider other systems where the conformal invariance is broken '4^. ^4ll |4^. [43l| . 
The model we presented in this paper is just a toy model. We need to investigate realistic models based on high 
energy physics which can accommodate our new mechanism. This would be an issue for future work. 
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